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Abstract.  Solar  wind  heat  flux  (HF)  electrons  are  valuable  as  tracers  of  the  interplanetary  magnetic  field  (IMF)  topology, 
distinguishing  positive  from  negative  solar  polarities  and  indicating  the  presence  of  magnetically  closed  CMEs  when  the  flows 
are  counterstreaming.  All  past  applications  of  heat  fluxes  to  determine  field  topologies  have  been  based  on  visual  inspection  of 
color  spectrograms  of  electron  pitch  angle  distributions  (PADs).  However,  HF  PADs  can  take  a  range  of  shapes  and  amplitudes, 
which  challenges  the  visual  analysis.  We  now  take  a  quantitative  approach  to  HF  analysis  by  parameterizing  the  HF  PADs 
of  the  UC  Berkeley  3DP  data  with  a  Fourier  harmonic  analysis.  We  have  calculated  the  harmonic  cosine  coefficients  Aq 
through  A4  for  a  five-year  period  of  the  Wind  3DP  data  set  with  a  10-min  time  resolution.  With  these  data  we  intend  to  derive 
quantitative  criteria  for  unidirectional  and  bi-directional  flows  and  other  possible  diagnostics  of  interplanetary  field  dynamics 
or  configurations.  Some  initial  considerations  and  results  of  the  3DP  parameterization  are  presented. 


1.  INTRODUCTION 

Solar- wind  heat  flux  electrons  have  provided  a  powerful 
tool  for  determining  the  solar  magnetic  polarities  of  the 
interplanetary  magnetic  field  (IMF)  [1,  2].  Those  elec¬ 
trons,  with  energies  E  >  80  eV,  stream  antisunward  par¬ 
allel  to  positive  polarity  field  lines  and  antiparallel  to 
negative  polarity  field  lines,  independently  of  the  local 
field  directions.  By  determining  whether  the  pitch  angle 
distribution  (PAD)  is  concentrated  at  0°  or  180°,  i.e.,  par¬ 
allel  or  antiparallel  to  the  field,  respectively,  we  can  es¬ 
tablish  the  solar  polarities  of  the  IMF.  It  is  possible  to  de¬ 
tect  fields  locally  turned  back  to  the  Sun  or,  with  bidirec¬ 
tional  electron  (BDE)  flows,  the  closed  fields  of  ICMEs 
(Figure  1). 

The  first  plots  of  HF  flow  directions  were  based  on 
a  solar-pointing  coordinate  system  and  required  a  visual 
comparison  of  the  observed  IMF  direction  with  the  HF 
flow  direction  to  estimate  the  PAD  (e.g.,  [3, 4]).  HF  plots 
from  current  spacecraft  are  done  in  a  coordinate  system 
fixed  on  the  IMF  direction,  which  immediately  yields  the 
PAD.  While  simplifying  the  analysis,  it  does  not  avoid 
the  considerable  variation  among  PADs  (e.g.,  [5,  6]) 
which  sometimes  leaves  the  net  HF  flow  direction  or 
the  presence  of  bidirectionality  [3]  in  doubt.  In  addition, 
the  PADs  are  usually  plotted  with  a  color  table  which 
may  or  may  not  be  normalized  for  the  total  number  of 
electrons  observed  in  the  PADs.  The  recent  reports  of 
HF  depletions  at  90°  PA  [6,  7]  suggest  another  possible 
HF  diagnostic  of  IMF  structure,  but  the  identification  and 
magnitude  of  these  features  is  compromised  when  we  are 


FIGURE  1.  Schematic  illustrations  of  open  fields  of  the  IMF 
locally  turned  back  to  the  Sun  (left)  and  a  closed  field  topology 
of  an  ICME  (right).  The  polarity  of  the  turned  back  fields  can 
be  inferred  from  the  primary  direction  of  the  HF  PAD,  and  the 
closed  fields  are  indicated  by  bidirectional  HF  PADs. 


limited  to  a  subjective  analysis  based  on  the  color-coded 
PAD  contour  plots. 


2.  ANALYSIS 

To  provide  quantitative  diagnostics  for  the  different 
types  of  PADs,  we  are  beginning  to  analyze  the  3DP 
[8]  HF  electron  PAD  data  in  terms  of  Fourier  cosine 
harmonics  using  the  following  least-squares  fit  to  the 


1,0 


0.5  - 

<t  0.0 - 

-0.5 

-1.nl .  .. 

0  90  180 

B  (deg) 


FIGURE  2.  Graphical  presentations  of  the  terms  of  the  fits 
to  equation  (1).  The  A2  and  A4  terms  are  important  for  bidirec¬ 
tionality  and  depletions,  respectively.  Dotted  lines  show  pro¬ 
files  of  negative  values. 


13-point  HF  PAD  intensities 

log  1(0)  =  Aq  +  Aicos(O)  -1-  A2Cos(20)  +  A3cos(30) 

+  A4Cos(40) 

where  6  is  the  PA,  which  ranges  from  0°  to  180®.  Be¬ 
cause  of  the  large  dynamic  range  of  the  HF  intensities, 
we  do  a  logarithmic  rather  than  a  linear  fit  to  the  inten¬ 
sity  1.  Plots  of  the  five  harmonic  terms  as  functions  of  9 
are  shown  in  Figure  2.  Note  that  a  bidirectional  flow  [3] 
should  have  a  significant  positive  A2  term  and  a  depletion 
[6,  7]  a  significant  negative  A4  term.  We  have  now  cal¬ 
culated  the  five  harmonic  components  for  all  the  WIND 
3DP  HF  data  from  1995  through  1999  in  10-min  inter¬ 
vals  for  electron  energies  of  125,  250,  and  500  eV. 

To  examine  the  signatures  of  the  various  An  we  se¬ 
lected  3  days  with  very  different  characteristic  PADs 
(Figure  3)  :  1  Jan  1995  has  a  poorly  defined  and  rela¬ 
tively  flat  PAD;  3  Aug  1996  has  a  well  defined  0®  PAD 
and  25  Jun  1998  has  a  BDE  PAD,  with  peaks  at  0®  and 
180®,  most  of  the  day. 


2.1.  A2  and  the  BDE  Diagnostic 

A2  is  the  key  parameter  for  detecting  BDE  flows. 
Figure  4  shows  the  plot  of  log  1(0)  for  various  ratios  of 
A2/A1  when  A2  and  Ai  >  0.  We  see  that  bidirectionality 
is  not  apparent  until  A2/A1  ^  1.  Selecting  an  arbitrary 
value  of  A2/A1  to  define  periods  of  bidirectionality  will 
be  a  critical  choice  for  our  anticipated  3DP  survey.  We 
note  that  Richardson  and  Reames  [9]  chose  A2/A1  >  0.8 
for  their  survey  of  bidirectional  energetic  (~  1  MeV)  ion 


FIGURE  3.  Three  days  with  very  different  3DP  260  eV 
electron  PADs:  top,  flat  PADs  of  1  Jan  1995;  middle,  directed 
PADS  of  3  Aug  1996;  bottom,  BDE  PADs  of  25  Jun  1998.  The 
PA  scales  extend  from  0®  on  the  bottom  to  180°  at  the  top. 


FIGURE  4.  Plot  of  log  1(0)  for  four  different  ratios  of 
|A2|/|Ai  I  with  other  A„  set  to  0. 


flows. 

Our  goal  here  is  to  get  a  general  impression  of  the 
parametric  variations  among  the  three  days  with  the 
very  different  PADs  as  determined  by  eye.  We  plotted 
log(|A2|)  against  log(|Ai|)  for  the  10-min  intervals  of 
each  of  the  three  selected  days  and  at  each  of  the  three 
electron  energies  used  in  the  harmonic  analysis.  Figure 
5  shows  the  results.  The  diagonal  lines  of  the  plots  en¬ 
able  us  to  compare  the  magnitudes  of  A2  with  A].  We 
see  first  that  the  plots  are  very  similar  for  each  of  the 
three  selected  electron  energies,  suggesting  that  the  par¬ 
ticular  energy  range  chosen  as  the  standard  for  analysis  is 
not  critical.  We  also  find  that  A2,  nearly  always  positive 
throughout  each  period,  exceeds  0.2  for  most  of  the  BDE 


FIGURE  5.  Plot  of  log(|A2|)  against  log(|Ai  |)  for  the  three 
different  electron  energies  and  for  the  three  different  PAD 
distributions  shown  in  Figure  3.  The  right  column  contains 
points  from  the  day  of  obvious  BDE  PADs.  Above  the  diagonal 
lines  IA2I  >  |Ai|. 


PAD  points,  but  is  rarely  that  large  in  the  flat  and  directed 
PADs.  In  addition,  while  A2  >  Ai  in  the  BDE  PADs  and 
A2  <  Ai  in  the  directed  PADs,  perhaps  as  expected,  we 
also  find  that  usually  A2  >  Ai  in  the  flat  PADs.  As  a  pos¬ 
sible  strategy  for  selecting  criteria  for  BDEs,  we  could 
plot  A2  versus  Aj  data  points  for  the  whole  1995-1999 
period,  looking  for  a  separate  population  of  points  in  the 
approximate  ranges  defined  by  A2  >  Ai  and  A2  >  0.2. 


2.2.  As  and  the  HF  Flow  Diagnostic 

The  most  important  goal  in  this  analysis  is  to  deter¬ 
mine  when  we  have  a  sufficiently  clearly  defined  HF  flow 
direction  to  decide  the  IMF  polarity.  From  Figure  2  we 
see  that  when  A3  and  Ai  have  the  same  sign  and  A3  > 
Ai ,  the  flow  direction  is  reversed  from  that  indicated  by 
Ai  alone.  In  Figure  6  we  show  plots  of  log(|A3|)  against 
log(|Ai|).  Points  above  the  diagonal  line  indicate  cases 
in  which  A3  is  large  enough  to  result  in  a  net  HF  flow 
reversed  from  that  of  Ai  if  both  terms  have  the  same 
sign.  A  substantial  number  of  points  in  the  flat  and  BDE 
PADs  suggest  such  reverse  flows.  However,  a  complica¬ 
tion  here  is  that  about  20-30%  of  the  A3  points  have  signs 
opposite  those  of  the  corresponding  Ai  values  and  there¬ 
fore  enhance  the  net  HF  flow.  In  addition,  large  values 
of  A2  could  act  to  diminish  the  net  HF  flows  determined 
from  the  Ai  and  A3  values  alone. 

The  ambiguities  involved  in  trying  to  sort  out  a  set 
of  relationships  among  Ai,  A2,  and  A3  suggest  that  a 
better  approach  may  be  simply  to  calculate  the  net  HF 


FIGURE  6.  Plots  of  log(|A3|)  against  logCjAiD  with  the 
same  format  as  in  Figure  5. 


over  all  electron  energies.  We  have  done  this  and 
show  the  plots  for  the  three  days  in  Figure  7.  The  red 
trace  is  parallel  to  the  magnetic  field  direction,  and 
the  purple  and  green  traces  show  the  Qe  in  directions 
orthogonal  to  the  IMF.  The  magnitudes  of  the  purple 
plots,  calculated  along  the  axis  perpendicular  to  the  plane 
of  the  IMF  direction  and  the  Sun,  are  comparable  to  the 
uncertainties  of  the  calculated  along  the  IMF  and 
can  be  used  as  filters  for  valid  HF  flow  directions.  In 
particular,  we  see  that  for  much  of  the  1  Jan  1995  flat 
PAD  plot  the  magnitudes  of  many  of  the  perpendicular 
HF  points  are  comparable  to  the  matching  points  of  the 
parallel  points,  suggesting  poorly  defined  HF  directions. 
On  the  other  hand,  the  parallel  HF  is  much  greater  than 
the  perpendicular  HFs  for  all  of  the  directed  PAD  plot  of 
3  August  1996.  Note  that  the  HF  calculation  used  here 
does  not  produce  the  systematic  20°  deviations  between 
the  Qe  and  B  directions  found  by  Salem  et  al.  [10]. 


2.3.  A4  and  Depletions 

The  inclusion  of  the  A4  harmonic  shown  in  Figure  2 
enables  us  to  look  for  the  cases  of  depletions  [6,  7],  the 
significant  decreases  symmetrically  centered  on  90°  PA. 
Figure  8  shows  plots  of  log(A4)  against  log( A2)  with  the 
positive  and  negative  values  of  log(A4)  plotted  on  sep¬ 
arate  panels.  Recall  that  A2  is  nearly  always  positive  in 
our  selected  data  sets,  producing  a  broad  PAD  depres¬ 
sion  at  90°.  A  negative  A4  term  is  consistent  with  the 
further  narrow,  symmetric  depression  at  90°  that  charac¬ 
terizes  depletions.  In  most  cases  of  Figure  8  IA4I  <  IA2I, 
and  A4  is  negative  (bottom  panels).  While  these  limited 


FIGURE  7.  Pairs  of  3DP  260  eV  PADs  (top)  and  Q,  (bot¬ 
tom)  for  the  same  three  days  of  Figure  3.  In  the  plots  the 
red  trace  shows  Q,  along  the  IMF  direction,  the  purple  along 
a  perpendicular  to  the  plane  of  the  IMF  direction  and  the  Sun 
and  the  green  along  the  third  orthogonal  axis. 


results  suggest  that  depletions  may  be  defined  by  the  cri¬ 
teria  that  A4  <  0  and  IA4I  >  0.5  x  IA2I,  a  much  larger 
survey  of  the  3DP  data  parameters  will  be  required  to 
justify  those  criteria. 


3.  DISCUSSION 

The  five-year  data  base  of  the  3DP  HF  electron  data  has 
been  analyzed  in  terms  of  the  first  four  Fourier  cosine 
coefficients  over  10-min  averages.  We  intend  to  develop 
quantitative  criteria  for  the  A„  coefficients  that  will  en¬ 
able  us  to  determine  IMF  solar  polarities  and  to  find  peri¬ 
ods  of  BDEs  and  depletions.  In  this  initial  work  we  have 
examined  characteristic  parametric  tradeoffs  among  the 
A„  for  three  selected  days  with  very  different  kinds  of 
PADS. 
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FIGURE  8.  Plots  of  log(A4)  versus  Iog(A2)  for  the  three 
days  of  the  study.  Here  we  distinguish  between  the  positive 
(upper  panels)  and  negative  (lower  panels)  signs  of  A4.  The 
predominately  negative  values  of  A4  are  expected  when  PAD 
depletions  [6, 7]  occur. 
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